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Abstract 

The meson spectra provide insight into the particle production mechanism and interaction in the 
hadronic and quark gluon plasma (QGP) phases. The detailed study of systematics of meson spec- 
tra is important also because it acts as ingredient for estimating the hadronic decay backgrounds 
in the photon, single lepton and dilepton spectra which are the penetrating probes of quark gluon 
plasma. In this work, we parameterize experimentally measured pion spectra and then obtain the 
spectra of other light mesons using a property known as mx scaling. The mj- scaled spectra for 
each meson is compared with experimental data for p+p, d+Au and Au+Au systems at ^/s]\^N 
= 200 GeV. The agreement of the m-r scaled and experimental data shapes are excellent in most 
cases and their fitted relative normalization gives ratio of meson to pion mx spectra. These ratios 
are useful to obtain the hadronic decay contribution in photonic and leptonic channels but also 
point to the quantitative changes in the dynamics of the heavy ion collision over p+p collisions. 
It is shown that, the particles with charm contents behave differently as compared to pions in 
d+Au systems and particles either with strange or charm contents behave differently from pions 
in Au+Au systems. For Au+Au system, three centrality classes have been studied which reveal 
that for the particles like kaon and (j), peripheral collision data is better reproduced as compared 
to central and their relative ratios with pions also increase as the collisions become more central. 
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I. INTRODUCTION 



Heavy ion collisions at relativistic energies are performed to create and study dense 
and/or hot matter in the laboratory. In Au+Au collisions at A/i^M_= 200 GeV at RHIC, 
many signals point to the formation of Quark Gluon Plasma (QGP) [ij. At RHIC, currently 
detailed properties of the strongly interacting matter are under investigation using a variety 
of observables. To isolate phenomena related to the dense and hot medium created in such 
collisions and to understand cold nuclear matter effects, it is also important to measure 
particle production in smaller collision systems like p+p and d+A. Measurements of trans- 
verse momentum spectra for particles emerging from p+p collisions are used as a baseline 
to which similar measurements from heavy ion collisions are compared. In addition, several 
observations from p+p collisions, such as the pt spectra with particle mass are interesting in 
their own right. The nuclear modification factor Raa for several identified hadrons at high 
transfer momentum serve as a probe for jet quenching [2|. It has been observed that there is 



strong suppression of hadrons in Au+Au collisions while no suppression is observed in d+Au 
collisions j3|. The photons Raa remains fiat both for d+Au as well as for Au+Au collisions 



even at high pt Being electromagnetic the photons escape unaffected and their spectra 
after subtracting the hadronic decay contribution reflect the properties of the medium such 
as temperature. The dielectron invariant mass measured by PHENIX collaboration shows 
enhancement over the cocktail from hadronic contribution in low mass region and is also 
associated .o the thennal radiation fro. QGP fl- The single electrons eo,n,ng fto,n sen,! 
electronic decays of charm and beauty quarks are important hard probes of the properties 
of matter produced in heavy ion collisions Q]. In order to obtain the electron spectra from 
charm decays one needs to subtract the electron contribution from other meson decays Q]- 
Getting the cocktail of single electrons, dielectrons and photon coming from the decays of 
all mesonic sources is crucial in any of the above analysis. 

WA80 collaboration found that the spectral shapes of vr and rj mesons are identical when 
plotted as a function of ttit js]. This property is known as ttit scaling and has been extremely 
useful to obtain the unknown meson spectra. 

In this work, we used itlt scaling to obtain all mesonic spectra from a given meson 
spectrum. We parameterize pion spectra first and then we obtain the spectra of other light 
mesons using tjit scaling. The relative normalization of the tjit scaled spectra is then fitted 



to the experimental data for all mesons in p+p, d+Au and Au+Au collisions. Both the 
magnitudes and shapes of the uit scaled and experimental data are studied for all these 
systems for many particles namely K, r;, 0, J/ifj and uj mesons. 



II. FIT PROCEDURE USING nix SCALING 

In this section, we describe the fitting procedure using rriT scaling, but before that we give 
a brief theoretical background of the fit function used in our analysis. Hagedorn proposed 



the following imperical formula jo] to describe the data of invariant cross section of hadrons 
as a function of pt over a wide range (0.3-10 GeV/c). Hagedorn described this as "inspired 
by QCD". 



j.d'N _ A 

Here A, po and n are fit parameters. The two limiting cases of this formula are as follows: 



— — ~ exp , for pj- — )■ (2) 

(1 + \ Po J ^ ^ 

~ ( — I , for Pt ^ oo. (3) 
\PtJ 

At low transverse momenta it assumes an exponential form and at large transverse momenta 
it becomes a power law arises from " QCD inspired" quark interchange model [l^ as: 

Even though the inclusive cross section is dominated by low pt particles, the hardening of 
the Pt distribution with increasing center of mass energy implies an increase of the transverse 
momentum. To better describe both low as well as high pt range, UAl collaboration [ll| 
used a hybrid form as follows: 



d'N 

E = B expi—brriT), for pt < Pr (5) 
dp'^ 

= (l^PT)n ^ for PT>P*T. (6) 
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3ENIX collaboration obtained a single form referred 
Q. The n,odmca.io„ . .o bet.e. d.cHbe the 



n. inis single lormuia nas been used to successiully describe 
in p+p [5|, Il2|-|l4| , d+Au [ij] as well as in Au+Au fl, Q, Q 



with free parameter. The PI 

as modified Hagedorn formula 
spectrum for wider pt range, in particular at high px where the spectrum behaves close 
to a simple power law function. This single formula has been used to successfully describe 
the hadron spectra measured 



collisions at different energies and is given by 

. ^ (,) 

dp' [exp(-apT - IvV) + 

which is close to a Hagedorn function at low pt and satisfies the requirement that the 
function is a pure power law at high pf . The term quadratic in p^ in the exponential term 
in the denominator is more important for the Au+Au collisions. 

In our analysis, first we parameterize experimentally measured pion spectra. In the 
fitting procedure, both neutral as well as charged pion data are incorporated. This is based 
on the assumption that the neutral pion spectrum is the same as the average charged pion 
spectrum. The fit function used here is the same modified Hagedorn distribution used by 
PHENIX collaboration but we replace pt by niT written as 



d^N A 
E- 



dp^ 



exp(— am^ — brnj.) + ^ 



UiJp'^ + ml), (8) 



where A, a, b, po and n are the fit parameters. The pion spectra measured in p+p, d+Au 
and Au+Au systems are fitted using this distribution and the parameters obtained are given 
in table [B For Au+Au system, data corresponding to three centrality classes namely 0-20 
%, 20-60 % and 60-92 % has been analyzed. The fit parameters for the centrality data 
for pions are given in table [Tll Comparing the values of b for different systems shows that 
term quadratic in rriT in the exponential term in the denominator is more important for 
the Au+Au collisions. We notice that the power n for all systems is close to 8 providing 
qualitative theoretical support for the quark interchange model. 

The light mesons which contribute sizeably to any measured electron and/or photons via 
their decay are pions, rj, p, u, 0, rjf etc. The r] meson contributes a sizeable fraction of decay 



electrons, in particular at high px- Using ttit scaling, we obtain the spectra of these light 
neutral mesons using pion fit function as 



(9) 



where rrih is the rest mass of the corresponding hadron or meson. The factor S is the 
relative normalization of the meson m-p spectrum to the pion m-p spectrum which we obtain 
by fitting the experimentally measured meson spectrum. 

The factor S should be close to, but not identical to, some standard meson/vr^ ratios in 



the literature 



15| that have been obtained as averaged over px intervals: 



r^/7r0=0.48 ± 0.03 [l6|; 



p/7rO=1.0 ± 0.3, predicted by [17|; 



w/7r°=0.9 ± 0.06 [isj; 

r]f/7T^=0A0 ± 0.12, predicted by PYTHIA |1 
0/7rO= 0.25 ± 0.08 |l9|; 



These values are obtained mostly by p+p measurements and models for p+p collisions, 
but also are widely used in cocktail calculations for the d+Au system. For Au+Au system 
the ratios are given in Ref. 20|. All the meson data used in this analysis along with mode of 
measurement and px ranges with their references are listed in table IIIII for p+p and in IIVI 
for d+Au and Au+Au systems. The particles like pions and kaons are measured by time 
of flight (TOF). The errors on the data are quadratic sums of statistical and uncorrelated 
systematic errors wherever available. All the data used in the present work is for same 
rapidity ( \y\ < 0.35 ) and from PNENIX experiment at RHIC. 



III. RESULTS 



Figure ((T^) shows the invariant yields of neutral 21| and charged pions 22| as a function 
of rriT measured in p+p collision at ^/snn = 200 GeV fitted with the Hagedorn function. 
The figure ([lb) shows the ratio of data to the fit. Figure ([2]) shows the invariant yields of 
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24| and J/ ip |25|, |26| as a function of measured in p+p system. 



The solid line is obtained using tjit scaling; the relative normalization has been used to fit 



TABLE I. The parameters of the Hagedorn distribution obtained by fitting pion spectra measured 
in p+p, d+Au and Au+Au colhsions at ^/snn = 200 GeV. 



Pa T a TTl Pf.PTS 


Colhsion systems 


p+p 


d+Au 


Au+Au 


A (GeV/c)-2 


11.11 ± 0.37 


52.30 ± 1.63 


822.30 ± 15.90 


a (GeV/c)-i 


0.32 ± 0.03 


0.24 ± 0.01 


0.420 ± 0.006 


b (GeV/c)-2 


0.024 ± 0.011 


0.11 ± 0.01 


0.215 ± 0.006 


Po (GeV/c) 


0.72 ± 0.03 


0.77 ± 0.02 


0.697 ± 0.003 


n 


8.42 ± 0.12 


8.46 ± 0.07 


8.35 ± 0.01 



TABLE II. The parameters of the Hagedorn distribution obtained by fitting pion spectra measured 
in Au+Au collisions of different centralities at y^s/vTV = 200 GeV. 



Parameters 


Au+Au cohision centrahty 


0-20 % 


20-60 % 


60-92 % 


A (GeV/c)-2 


1883.23 ± 31.20 


744.77 ± 9.11 


132.99 ± 2.21 


a (GeV/c)-i 


0.442 ± 0.006 


0.385 ± 0.004 


0.273 ± 0.005 


b (GeV/c)-2 


0.242 ± 0.006 


0.192 ± 0.004 


0.154 ± 0.006 


PO (GeV/c) 


0.708 ± 0.003 


0.694 ± 0.003 


0.653 ± 0.006 


n 


8.40 ± 0.01 


8.26 ± 0.01 


8.14 ± 0.03 



the measured spectra. Figure ([3]) shows the invariant yield of u meson 13|, l27|| as a function 
of rriT measured in p+p system along with the rriT scaled curve (solid line). 



Figure shows the invariant yields of neutral [28| and charged pions [22| as a function 
of rriT measured in d+Au at a/J/vtv = 200 GeV fitted with the Hagedorn function. The 
figu re (Bb ) shows the ratio of data to fit. Figure ([5]) shows the invariant yields of measured 



K± |22|], [16!,|28l, |24| and J/tP [29] function of rriT measured in d+Au system. The 
solid line is obtained using rriT scaling; the relative normalization has been used to fit the 
measured spectra. Figure (EI) shows the invariant yield of u meson as a function of rriT 
measured in d+Au system along with scaled curve. 



Figure ([7^) shows the invariant yields of neutral 3l| and charged pions 3^ as a function 
of niT measured in Au+Au at ^/sNN = 200 GeV fitted with the Hagedorn function. The 



TABLE III. Particles with their measured decay channels and pt range for different p+p collisions 
with references. The data is for central rapidity region \y\ < 0.35. 



Particle 


Mode 


Pt range 


Reference 


vrO 


/ / 


6-18 9 GeV 


[211 


7r± 


TOF 


3-2 6 GeV 


[221 


K± 


TOF 


0.4-1.8 GeV 


[221 


71 
'1 


1 I 


2 7-11 GeV 


[161 






5 5-21 GeV 


[231 


T 




0-3.5 GeV 


[13.1 






1-7.0 GeV 


[241 


CO 




0.1-3.5 GeV 


[13j 






2.5-9.3 GeV 


[22 






2.25-13 GeV 


[13j 






0.1-8.5 GeV 


[25] 






0.1-8.5 GeV 


[26] 



figure 



K 



± 



Uh) shows the ratio of data to fit. Figure (IHl) shows the invariant yields of measured 



32|, 7] 23|, |30|, (j) and J/ip [33*] as a function of rriT in Au+Au system. The solid line 
is obtained using rriT scaling; the relative normalization has been used to fit the measured 
spectra. 

n 

Figure ([9]) shows the invariant yield of u [3j] meson as a function of rriT measured in 
Au+Au system along with mj- scaled curve. 

The table IVl shows the normalization factor of meson to pion ttlt spectra (meson/7r°) 
obtained by fitting the m-r scaled spectra with the measured spectra for different collision 
systems. In case of p+p they are in agreement with the ratios available in the literature. 
One can observe that the shapes of the derived uit scaled spectra very well reproduce the 
measured spectra. 

In case of d+Au, the fitted meson to pion ratios are in agreement with those in case of 
p+p for all mesons except for the The J/ip to pion ratio is small in d+Au case because 
the J/ip yields are suppressed but the pion yields remain unaffected. The shapes of the 
derived mx scaled spectra reproduce all measured spectra including J/ip. It shows that the 



TABLE IV. Particles with their measured decay channels and px range for d+Au and Au+Au 
systems with references. The data is for central rapidity region \y\ < 0.35. 



Particle 


Mode 


range 


Reference 


d+Au collision 


/I 


77 


19 17 PpV 
1..^-! 1 .u vje V 






TOF 


0.3-2.6 GeV 


[22] 




TOF 


0.45-1.8 GeV 


[22j 


V 


77 


1.2-4.75 GeV 


[16] 






2.25-11.0 GeV 


[28] 




K+K- 


1.1-7 GeV 


m 


w 




3-9 GeV 


m 




e+e^ 


0.5-4.5 GeV 


[29] 


Au +Au collision 




77 


1.2-19.0 GeV 


m 




TOF 


0.25-2.95 GeV 


[32] 




TOF 


0.45-1.95 GeV 


[32j 


V 


77 


2.25-9.5 GeV 


[30] 






5.5-21 GeV 


[23] 




K+K- 


1.1-7 GeV 




OJ 




4.5-8.5 GeV 


[34] 






0.25-5 GeV 


[33] 



^/ip yields are suppressed uniformly in all the pt range considered here. 

In case of Au+Au, the fitted 7] to pion ratio is similar to that in case of pp. It is because 
7] and pions are suppressed by similar amount. Also the shape of the curve matches well 
with the data. The same is true with uj. It means that one can obtain the ratios of the 
and u rriT spectra with pion m-r spectra for p+p collisions and then can use them for the 
case of d+Au and Au+Au systems. The kaon to pion ratio is larger than that in p+p and 
the rriT scaled curve does not reproduce the shape of the measured spectra that well. The 
to pion ratio is larger than that in p+p because the is not suppressed as much as pions. 
One can observe in case of that the shape of measured spectra is not well reproduced 
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TABLE V. The relative normalization (meson/vr") obtained by fitting the my scaled spectra with 
the measured spectra for different collision systems. 





Collision systems 


P+P 


d+Au 


Au+Au (MB) 




0.422 ± 0.003 


0.439 ± 0.002 


0.530 ± 0.002 


r//7r0 


0.497 ± 0.015 


0.460 ± 0.017 


0.525 ± 0.040 




0.233 ± 0.011 


0.215 ± 0.008 


0.348 ± 0.017 




0.903 ± 0.021 


0.964 ± 0.10 


0.845 ± 0.145 




0.054 ± 0.002 


0.0021 ± 0.0002 


0.0037 ± 0.0003 



TABLE VI. The relative normalization (meson/vr") obtained by fitting the m-r scaled spectra with 
the measured spectra for different Au+Au collision of different centralities. 



particle ratio S 


Au+Au collision 


MB 


0-20% 


20-60% 


60-92% 


K/vrO 


0.530 ± 0.002 


0.504 ± 0.001 


0.494 ± 0.001 


0.486 ± 0.003 




0.525 ± 0.040 


0.542 ± 0.017 


0.580 ± 0.014 


0.548 ± 0.024 




0.348 ± 0.017 


0.401 ± 0.010 


0.391 ± 0.007 


0.302 ± 0.014 




MB 


0-20% 


20-40% 


40-92% 


0.0037 ± 0.0003 


0.0032 ± 0.0005 


0.0043 ± 0.0005 


0.0033 ± 0.0006 



in intermediate rriT region of — m ~ 2 to 4 GeV.. This shows that the particles with 
strangeness contents behave differently in Au+Au case. 

The J/^ to pion ratio is small in Au+Au case but even smaller in d+Au. This is quite 
interesting and is because the yield of pions is suppressed only in Au+Au but the yield 
of i/ip is suppressed in both systems. Interestingly, the shape of the derived rriT scaled 
spectra reproduce the measured i/ip spectrum. It shows that the i/ip yields are suppressed 
uniformly in all the mx range considered here. 

We have also considered three centrality classes for Au+Au system; 0-20 % (most central), 
20-60 % (intermediate) and 60-92 % (peripheral). In case of J/?/) we have data available for 
the centralities 0-20 %, 20-40 % and 40-92 %. The class 40-92 % in this case can be called 
as semiperipheral. Also we do not have good centrality data for oj. The table IVl shows the 
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normalization factor of meson to pion ttit spectra (meson/vr") obtained by fitting the tut 
scaled spectra with the measured spectra for different centralities. 

Figure [10] shows the invariant yield of measured neutral 3l| and charged pions 32| as a 
function of in Au+Au system for different centralities. The solid lines are the Hagedorn 
function fits, the parameters of which are given in Table [Tll Figure [TT] shows the invariant 
yield of measured 32|] as a function of rriT in Au+Au system for different centralities. 
The solid lines are obtained using ttlt scaling from the corresponding centrality data of 
pions. The relative normalization has been used to fit the measured spectra. Figure [T2] 
shows the invariant yield of measured r] (open symbols 30|, solid symbols 23|) as a function 
of rriT in Au+Au system for different centralities. Figure [13] shows the invariant yield of 
measured as a function of rriT in Au+Au system for different centralities. Figure [H] 
shows the invariant yield of measured J/ip 33| as a function of in Au+Au system for 
different centralities. 

In case of r] the data of all centralities are very well reproduced by m-p scaled pion data. 
The ratios obtained also remain approximately same. In case of kaons, the shapes of the 
peripheral data are reproduced by rriT scaled pion data very well but for the most central 
collisions the disagreement between the two is quite evident. Also the fitted ratios tend to 
increase as we move from peripheral to central collisions. The same seems true for where 
peripheral data is better reproduced as compared to central data and the ratios also increase 
as the collision becomes more central. In case of J/ip all three centralities seem same both 
by shape and relative ratio to the pions. At least from the present data they can not be 
distinguished. 



IV. CONCLUSION 

In summary, we parameterize experimentally measured pion spectra and then obtain the 
spectra of other light mesons using a property known as m-r scaling. The m-r scaled spectra 
for each meson is compared with experimental data for p + p, d + Au and Au+Au systems 
at ^ysNN = 200 GeV. Their fitted relative normalization gives meson to pion ratio. These 
ratios would be useful to obtain the hadronic decay contribution in photonic and leptonic 
channels. The shape of the derived mj- scaled spectra very well reproduce the measured 
spectra and fitted meson to pions are in agreement with the ratios available in literature 
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for p+p system. We show that the ratios of the t] and u ttlt spectra with pion rriT spectra 
obtained for p+p coUisions can be used for the case of d+Au and Au+Au systems. It is 
observed that the particles with charm contents behave differently from pions for both the 
d+Au and Au+Au systems. The particles with strange contents behave similar as pions 
in d+Au system but behave differently from pions in Au+Au systems. The more detailed 
centrality analysis of Au+Au collision reveal that in case of r] the data of all centralities are 
very well reproduced by itlt scaled pion data. The ratios obtained also remain approximately 
same. For particles, kaon and 0, peripheral data is better reproduced as compared to central 
data and the ratios also increase as the collision becomes more central. In case of J/^/', all 
three centralities seems same both by shape and relative ratios to the pions within the errors. 
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FIG. 1. (Color online) a) the invariant yields of neutral [21!] and charged pions 
of rriT measured in p+p collision at 200 GeV fitted with the Hagedorn function, b) the ratio of 
data to the fit. 
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p+p Minimum Bias @ 200GeV 
□ <!>-> e*e' 
• (|)->K*K" 
— m^ scaied curve 

Fitted norm. = 0.233 ± 0.011 




10" 



o 10' 



p+p Minimum Bias @ 200GeV 
□ J/T ^ e*e" 
• J/'P->eV 
— m^ scaled curve 

Fitted norm. = 0.0543 ± 0.0016 



2 3 4 5 
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2 3 4 5 
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23|), 



FIG. 2. (Color online) The invariant yield of 22], (o pen squares [l^, solid circles 
(open squares [3], solid circles 24 1) and J (open squares 25], solid circles [2^) as a function of 
rriT measured in p+p system. The solid line is obtained using rriT scaling; the relative normalization 
has been used to fit the measured spectra. 



14 




d+Au Minimum Bias @ 200GeV 

A (7t* + 7t")/2.0 




6 8 10 12 14 16 18 
- m [GeV/c^] 



2 
1.8 
1.6 
1.4 

1.2 

0.8 r 
0.6 
0.4 
0.2 
0, 



(b) 



n 



d+Au IVIinimum Bias @ 200GeV 
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FIG. 4. (Color online) a) the invariant yields of neutral [28] and charged pions 
of rriT measured in d+Au at 200 GeV fitted with the Hagedorn function, b) the ratio of data to 
fit. 
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?IG. 5. (Color online) The invariant yield of measured [22], rj (open squares 16 1, solid circles 
28l|). (p 241 and J/ip [29] as a function of mr measured in d+Au system. The solid line is obtained 



using rriT scaling; the relative normalization has been used to fit the measured spectra. 
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FIG. 7. (Color online) a) the invariant yields of neutral [31'] and charged pions 32] as a function 
of rriT measured in Au+Au at 200 GeV fitted with the Hagedorn function, b) the ratio of data to 
fit. 
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FIG. 9. (Color online) The invariant yield of u 
system along with nriT scaled curve. 
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FIG. 10. (Color online) The invariant yield of measured neutral [31] and charged pions [32] as a 
function of mx in Au+Au system for different centralities. The solid lines are the Hagedorn fit 
function. 
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FIG. 11. (Color online) The invariant yield of measured clS cl function of m-p in Au+Au 

system for different centralities. The solid lines are obtained using rriT scaling; the relative nor- 
malization has been used to fit the measured spectra. 
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FIG. 12. (Color online) The invariant yield of measured r/ (open squares 
a function of rriT in Au+Au system for different centralities. The solid lines are obtained using my 
scaling; the relative normalization has been used to fit the measured spectra. 
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FIG. 13. (Color online) The invariant yield of measured (j) 241] as a function of mx in Au+Au system 
for different centralities. The solid lines are obtained using my scaling; the relative normahzation 
has been used to fit the measured spectra. 



0.5 



Au+Au @ 200GeV = 
J/y -> e*e' 
Minimum Bias x 10 
0-20% centrality x 2 E 
20-40% centrality I 
40-92% centrality 
scaled curves E 




1 1.5 2 
mj-m [GeV/c^] 



2.5 



function of my in Au+Au 



FIG. 14. (Color online) The invariant yield of measured J/ip 
system for different centralities. The solid lines are obtained using m-p scaling; the relative nor- 
malization has been used to fit the measured spectra. 
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